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In this study, we have performed a detailed investigation of the electronic properties 
of a core/shell/well/shell mult i- layered spherical quantum dot, such as energy eigen- 
values, wave functions, electron probability distribution and binding energies. The 
energy eigenvalues and their wave functions of the considered structure have been 
calculated for cases with and without an on-center impurity. For this purpose, the 
Schrodinger equation has been numerically solved by using the shooting method in 
the effective mass approximation for a finite confining potential. The electronic prop- 
erties have been examined for different core radii, barrier thickness and well widths 
in a certain potential. The results have been analyzed in detail as a function of the 
layer thickness and their physical reasons have been interpreted. It has found that 
the electronic properties are strongly depending on the layer thickness. 
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I. INTRODUCTION 



The recent developments in the manufacturing and growth technology have caused the 
fabrication of the low dimensional systems which are called nanostructures. These structures 
are different from the macroscopic structures and they present rather interesting properties.- 
Because of these properties, the nanostructures have a wide area of the research and interest 
for the scientists.- The electron motion can be confined in all three dimensions with advanced 
manufacturing technology and so, these improvements have made possible the production 
of quantum dots (QDs). These confinements have led to sharp energy levels and density 
of states as which were found in atoms. Because of that, QDs are also called artificial or 
super atom.« Outstanding properties of QDs are changes in energy spectrum by doing some 
modifications in systems which have different geometry and formed via growth technology. 
The investigations claim that the electronic properties of QDs are seriously connected to 
their size, shape and dimensions.- Because of the possibility to control of the sizes, shapes, 
energy levels and electron number of QDs, there has been considerable amount of theoretical 
and experimental studies on these structures. 

As is well known, an impurity atom in a semiconductor is a different one from the semicon- 
ductor's atom and it has an additional electron more than that required to make a chemical 
bonds with neighboring atoms. This more electron can be easily donating to the host ma- 
terial. Such impurity atoms are known as donors. If a donor impurity has an additional 
one electron and is unionized, this kind of impurity is known as hydrogenic one because of 
that it is very similar to a hydrogen atom. Thus, most of properties of this impurity are 
overlap with that of a hydrogen atom. An example of a typical hydrogenic donor impurity 
in a GaAs crystal is silicon.- More detail related with a hydrogenic donor impurity can be 
found in Ref. 

The existence of impurities causes some changes in the performance of quantum devices 
and their electronic, optical and transport properties because of the Coulomb interaction 
between the electron and the impurity atom. Although it does not provide fully accurate 
descriptions about the QDs, the study on hydrogenic impurity problem in semiconductor 
QDs is very useful model, for understanding the optical and the electronic properties of 
these structures. Therefore, this model is used in a number of studies, since many physical 
properties of an impurity can be explained in these kinds of heterostructures, even though 

2 



it is simple and a somewhat rough estimate.— _ — The confinement in all three dimensions 
in QDs causes a reduction of the distance between the electron and the impurity leading 
an increase in the Coulomb interaction. The hydrogenic impurity states in the QDs have 
been widely studied by many authors due to their importance for potential applications 
in optoelectronic and photonic devices. Zhu at al.— and Porras-Montenegro and Perez- 
Merchancano^ have carried out pioneering studies about confinement effect on hydrogenic 
impurity states in QDs. In these studies, binding energy has been probed as a function of 
QD's size and impurity location. In addition to these studies, the hydrogenic donor binding 
energy has been surveyed as theoretically by many authors for QDs which have different 
geometries such as spherical, cylindrical, pyramidal and rectangular.— r — Riberio and LatgeSi 
worked on a cubic quantum heterostructures and they have showed that the binding energy 
values are close to those of cubic and spherical QDs which have near volume. Bose and 
Sarkar— have used perturbation and variation technique to obtain hydrogenic donor binding 
energies in a spherical QD for an infinite parabolic confinement. The study of the impurity 
effect on the binding energy has made possible the fabrication of a low threshold laser diode. 23 
Jayam 2 ^ has investigated the effect of the electric field on a hydrogenic donor binding energy 
and same problem has been studied by Peter— for the the effect of a hydrostatic pressure. 
Hydrogenic donor binding energy for excited states in a spherical QD has been analyzed by 
Sadeghi.^ 

Recent improvements in growth technology have resulted in manufacturing of multi- 
layered spherical quantum dots (MSQD) and so, MSQD has been studied theoretically 
by many researchers.— The binding energy of a hydrogenic impurity in a GaAs-AlGaAs 
MSQD has been calculated as a function of barrier and the inner dot thickness for differ- 
ent potentials. — For the same structure, Boz et al.— explored effects of the geometry on 
hydrogenic impurity energy states by using a fourth order Runge-Kutta method. 

In this paper, we investigate the electronic properties of a core/shell/ well/shell QD struc- 
ture for both ground (Is) and excited (lp) states. We compare these properties for cases 
with and without the impurity and for different layer thickness. At the same time, the effect 
of the variations in the geometry of the structure on the donor binding energy is studied. 
All results are presented for both cases with and without the impurity. 

The rest of the paper is organized as follows: In the next section, we describe the model 
and calculations briefly. Results and discussion are presented as detailed in section III. In 
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the last section, a short conclusion is given. 



II. MODEL AND CALCULATIONS 

In this study, we have considered a spherical core/shell/well/shell QD containing a single 
electron. This structure comprises of two GaAs dots one embedded in the other. Also, in 
this structure, each of QDs is isolated by means of AlGaAs which provides a finite confining 
potential for the electron. Out of this structure has been coated by AlGaAs again as seen 
in FigJTJ In our system, spherical QDs have been interacted through potential barrier 
which allows electrons to tunnel between the core and well regions. The shell thickness 
is T s = i?2 — R\ and the well width is T w = R% — R 2 . The effective mass approximation 
and BenDaniel-Duke boundary conditions are used in the calculations. Here, the impurity 
position is taken as center of the structure. In the effective mass approximation, the single- 
particle Schrodinger equation is given as 

f-^V r ) - — + + V b ] R n , (r) = E n/ R n , (r) . (1) 

Here, h is reduced Planck constant, m*(r) is the position-dependent electron's effective 
mass, Z is the charge of the impurity, i is the angular momentum quantum number, Vb is 
the finite confining potential, E n g is the electron energy eigenvalue for a specified principle 
quantum number (n) and angular momentum quantum number (£), R n ^(r) is the radial 
wave function of the electron. It should be noted that Z = and Z = 1 correspond to cases 
without and with a hydrogenic donor impurity, respectively. The mathematical expression 
of the confining potential is 

0, r <R l and R 2 < r < R 3 

V h =\ . (2) 

228 meV, R x < r < R 2 and r > R 3 

Here, 228 meV is the conduction band offset in GaAs/AlGaAs junction.—^ 

In order to determine the single particle energy levels and corresponding wave functions, 
Eq.([TT) is solved full numerically. Shooting method is used for the determination of the 
single particle energies and corresponding wave functions. As is well known, this technique 
converts an eigenvalue problem to an initial value problem. For this purpose, Hamiltonian 



FIG. 1. Schematic representation of a core/shell/well/shell QD and its potential profile. 

operator is discretized on a uniform radial mesh in ID using the finite differences, then 
Eq.flT]) can be reduced to an initial value equation by means of 
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where i is the index of mesh points, h is the distance between two mesh points and it is 

n 

chosen as 0.005. The details of this method can be found in Ref. [9|. 

The neutral donor binding energy depends on position of the impurity and geometry of 
the structure 30 and is determined as 



E b = E n> i — E U; e(D 



(4) 



where E n> £ and E n ^(D°) are the energy levels for cases without and with the impurity, 
respectively. 



III. RESULTS AND DISCUSSION 



The atomic units have been used throughout the calculations, where Planck constant 
h — 1, the electronic charge e = 1 and the bare electron mass itlq = 1. Effective Bohr radius 
is a = 0.529-^ — 100A and effective Rydberg energy is R v = 13.6^- ~ 5.25 meV— 
The material parameters have been taken as^-^ mG a A s =0.067m , m J 4;GaAs=0.088m , 
K Ga A s =l3.18 and KAiGaAs=^-^- Also the effective masses of electrons inside GaAs and 
AlGaAs are m* and m^, and the dielectric constants are K\ and K2, respectively. The 
position-dependent effective mass and the dielectric constant may be defined as follows 33 

5 




FIG. 2. Ground and excited state energy values as a function of core radius for constant shell 
thickness (T s = 0.2ao) and two different well widths (T w = 0.4 and 0.6ao) for Z = (top panel) 
and Z = 1 (bottom panel). 
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1, r <R l and R 2 < r < R 3 
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3, R l <r<R 2 and r > R 3 
1, r < i?i and R 2 < r < R 3 

^, Ri<r < R 2 and r > R 3 



(5) 



A. Effect of the layer thickness on the electronic properties 

Figure [2] shows the variation of the energy levels for Z = (top panel) and Z = 1 (bottom 
panel) as a depending upon Ri for both ground (Is) and excited (lp) states in a MSQD. 
In this calculation, the barrier thickness is taken as T s = R 2 — Ri = 0.2a and the results 
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FIG. 3. The electron density of ground and excited states for Z = (top panel) and Z = 1 (bottom 
panel). The thickness of the layers are specified on the figures. 

are comparatively given for two different well widths specified on the figures. As seen from 
the figures, when R\ increases, both ground and excited state energies exhibit decreasing 
tendency. This is an expected result because of the relation between energy and radius, 
E oc 1/R\. Also, we observe that ground and excited state energies have different values 
until certain Ri value for different well widths in both Z = and Z — 1 cases. After this 
Ri value, ground states energies close to each other and they have same values with further 
increasing of the core radius. At the same time, the figures show that these consequences 
are valid for also excited states. In Z = case, while ground state energies have same values 
at around R\ = l.Oao, these changes are at about R\ = lAao for excited state ones. In 
case of Z = 1 , this behaviors are observed at approximately R\ = 0.7ao and R\ = 1.2ao for 
ground and excited states, respectively. This is because, in small value of the core radius, 
probability of ground and excited state electrons in the well region is larger than that in the 
core region as seen from Fig. [3j When the core radius increases, the electron probability 
starts to increase in this region due to the rule of energy minimum tendency and relation 
between energy and core radius as mentioned above. This increase in the probability is more 
evident for case with Z = 1 due to attractive Coulomb potential of the impurity as seen 
in Fig. [3] (e). After the specified core radii, the probability densities of both ground and 
excited states with and without the impurity become maximum in the core region as seen 
in Fig. |3] (c) and (f). Because, the electron and energy levels do not feel the effect of the 
well region. As a result, core/shell/well/shell structure exhibits the properties of a single 
QD with a finite potential barrier. 
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In the system under consideration, AlGaAs layer embedded between two GaAs QDs 
creates a finite potential barrier and so, the electron tunneling occurs between the core and 
well regions. Thus, it can be expected that the variation of the barrier thickness affects the 
energy values of this system for cases with and without the impurity. In order to understand 
these variations, we must analyze Fig. H] which shows the variation of the energies of both 
ground and excited state energies with the shell thickness for Z = (top panel) and Z = 1 
(bottom panel). The figures are plotted for two different core radii (Ri = 0.8 and l.Oao) 
and same well widths (T w = 0.2ao). As seen from these figures, in all cases, while the energy 
values increase until a small fraction of the shell thickness, they remain constant values after 
those thickness. These tendencies can be explained as follows: When the barrier thickness 
is small, the wave functions of ground and excited states can tunnel between the core and 
well region for both Z = and Z — 1 cases as seen in Fig. |5j As a result of this, the 
electron which is in the core region (or in the well region) is affected by the well region (or 
core region) for both ground and excited states. This response is less evident for Z = 1 
case because of attractive Coulomb potential of the impurity. However, an increase in of the 
barrier thickness decreases possibility of the tunneling between core and well regions and so, 
the electron confined in the core region can not sense the effect of the well region. Because 
of that, the energies have same value in the case of further increasing of barrier thickness. 
We conclude that the barrier thickness is rather important until a critical value for tuning 
the electronic structure of a MSQD. 

Fig. [6] shows the variation of the ground and excited state energy values with the well 
width (T w ) for Z = (top panel) and Z = 1 (bottom panel). The figures are plotted 
for two different barrier thicknesses (T s = 0.2 and 0.4ao) and the core radius is taken as 
Ri = 0.6ao- As seen from the figure, the energy levels, especially Is states, remain constant 
until T w = OAclq for case without impurity, when T s = 0.4a . This treatment continues 
to approximately T w = 0.6ao for case with impurity because of the attractive Coulomb 
potential of the impurity. These changes are more smooth for T s = 0.2ao, because the 
electron is localize in the core region (or in the well region), it can leak easily to the well 
region (or to the core region). These treatments can be obviously observed in Fig. [7J In the 
both cases of Z = and Z — 1, the difference between ground and excited state energies is 
considerably large for the small well widths. All the energies decrease with further increase in 
the well widths. As is well known, when the well width increases, the quantum confinement 
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FIG. 4. Variation of ground and excited state energy values with shell thickness (T s ) for constant 
well width (T w = 0.2ao) and two different core radii = 0.8 and l.Oao) for Z = (top panel) 
and Z = 1 (bottom panel). 
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FIG. 5. The electron density of ground and excited states for Z = and two different core radii, 
R\ = 0.8ao (top panel) and R\ = l.Oao (bottom panel). The thickness of the layers are specified 
on the figures. 
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FIG. 6. Variation of ground and excited state energy values as a function of well width (T w ) for 
constant core radius = 0.6ao) and two different shell thickness (T s = 0.2 and 0.4ao) for Z = 
(top panel) and Z = 1 (bottom panel). 



effect decreases. Therefore, the consecutive energy levels become close to each other in both 
Z = and Z = 1 cases. In addition, as seen from Fig. the electron escapes from the 
core region and it is completely confined in the well region with increasing well widths. If 
the well region is much more extended, the electron feels less confinement effect on it and 
hence, the energy values of the electron go towards a free particle ones. Moreover, ground 
and excited states energies for Z = 1 case are smaller than those for Z = case. Also, the 
difference between ground and excited states energies in Z = 1 case is larger than those for 
Z = case. All of these situations are originated from the attractive Coulomb interaction 
between the electron and impurity. 
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FIG. 7. The electron density of ground and excited states for Z = and two different well width, 
T s = 0.2ao (top panel) and T s = 0.4ao (bottom panel). The thickness of the layers are specified on 
the figures. 
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FIG. 8. Variation of ground and first excited state impurity binding energies with core radius {R\) 
for T w = 0.6ao and two different shell thickness, T s = 0.2 and 0.4ao- 

B. Effect of the layer thickness on the donor binding energy 

Variation of ground and first excited state binding energies of the donor impurity as a 
function of the core radius {R\) is given in Fig. [8] for constant well width and two different 
shell thickness that are specified on the figure. As seen from the figure, ground state binding 
energy exhibits a decreasing tendency until almost R\ = 0.4ao value for T s = 0.2ao- This 
decreasing tendency proceeds until R\ = 0.9ao for excited state. This treatments can be 
explained as follows: In smaller values of core radius, the electron density in the well region 
is larger than that in the core of the structure for both Z = and Z = 1 cases. Hence, the 
attractive Coulomb potential of the impurity does not pull so much down the energy levels. 
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FIG. 9. Variation of the ground and excited states electron densities with r for constant shell 
thickness and well widths for Z = (top panel) and Z = 1 (bottom panel) cases. 

The expanding core region, however, is to become much more effective on the decreasing of 
energy levels. After certain Ri values, the binding energies start to increase rapidly until 
Ri = 0.9ao for Is state. In case of lp state, this increase, ongoing until R± = 1.2ao, is more 
smooth compared to that in ground state. This behavior indicates that the probability of 
electron tunneling to core region rises with increasing of the core radius. In this case, the 
impurity pulls the energy values more down, especially ground state energies as expected, and 
so, this circumstance causes an increase on the donor binding energy. With further increasing 
of the core radius, the binding energies decrease smoothly and go towards bulk values due 
to the vanishing of the confining effects on the electron. Although similar behaviors are also 
observed for T s = 0.4ao case, the binding energy is bigger and their increasing becomes too 
steep. As is well known, if the barrier thickness increases, the tunneling probability of an 
electron confined in the core region (the well region) to the well region (the core region) 
decreases. We can easily conclude that the shell thickness has an important effect on the 
binding energies. 

In order to make clear the impurity effects on the electron densities of ground and excited 
states, we plot Figj9] for different core radii and constant shell thickness and well widths. 
Figures M (a) and (d) show the electron densities in the same core radius for Z = and Z = 1 
cases, respectively. As seen from Figj9](a), although the ground state density is substantially 
in the well region, the excited state density is almost completely in the well layer. When we 
compare Fig. [9] (d) with Figj9] (a), we can see that there is almost no effect of the impurity 
on the excited state density. However, the ground state density is seriously affected from 
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FIG. 10. Ground and first excited states impurity binding energy of a core/shell/ well/shell QD 
as a function of the shell thickness (T s ) for T w = 0.2ao and two different core radii, R± = 0.8 and 
1.0a . 

the attractive coulomb potential of the impurity and as a result of this, the ground state 
density becomes localize substantially in the core region. We conclude that although the 
impurity has no serious effect on the excited state energies, it is much more effective on the 
ground states and it pulls the ground state energies more down. In Figj9](b) and (e), it can 
be seen that the ground state densities become localize in the core region with increasing 
core radius for cases with and without the impurity. But, the strength of this localization 
is a bit more for Z = 1 case, when it is compared to the Z = case. In addition, we can 
see that the impurity has somewhat effect on the excited states, even though it is lesser. In 
Figj9] (c) and (f), when the core radius R\ = 1.5a , both ground and excited state densities 
are almost completely localize in the core region and the core/shell/well/shell QD structure 
exhibits a single core/shell QD properties with further increasing of the core radius. 

In multi-layered spherical quantum dot, the change in the barrier thickness is effective 
on the binding energy till a critical value of the thickness due to tunneling, which depends 
strongly on the barrier thickness, between the core and well region. The variation of the 
binding energy with the shell thickness is plotted in FigJlO] for two different core radii 
= 0.8 and 1.0a ) and well width T w = 0.2ao- As seen from the figure, the binding 
energy of ground state increases until T s = 0.3ao and it reaches up a certain value with 
further increasing of the T s for Ri = 0.8ao- This increase is continuing up to T s = 0.5ao 
for Ri = l.Oao- This is that because the core region becomes large, the electron feels lesser 
the effect of the well region. In excited states binding energies, these increases are observed 
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FIG. 11. Variation of the impurity binding energy of ground and first excited states with the well 
width (T w ) for T s = 0.4ao and two different core radii, i?i = 0.6 and 0.8ao- 

until T s = 0.2a value for both R\ = 0.8 and l.Oao- This behavior in binding energy can 
be interpreted as follows: In small shell thickness and when the well width is narrow, the 
electron has been confined into the core region. In addition to this, because the barrier 
thickness is thin, the effect of the well has been felt by electron. Therefore, as well as 
the impurity, the well has an influence on the energy levels and plays an important role in 
reducing the binding energy. As a result, in this shell thickness, the binding energy is small. 
If the barrier thickness exceeds a certain critical value, the electron in the core no longer feel 
the effect of well region and, therefore, the impurity binding energy remains constant with 
increasing shell thickness. Similar process is valid for excited state binding energies. 

Figure [11] shows the variation of the donor binding energy with the well width for two 
different core radius (Ri = 0.6 and 0.8ao) and shell thickness T s = 0Aa . If we look at the 
figure, the donor binding energies remain almost constant values for both ground and excited 
states (except for case with R\ = 0.6ao) until certain well widths. After this certain values, 
these binding energies decrease sharply and approach to the same values. The physics of 
these behaviors can be explained as follows: When the well width is small, the electron 
becomes localize in the core region. Therefore, the impurity pulls more down the ground 
state energy and causes higher binding energies. With increasing well width, the electron 
in the core starts feeling the well region and while the probability distribution decreases in 
the core, that increases in the well. As a result of this, the binding energy decreases with 
increasing well width. If the width is increased further, the electron in the well does not feel 
the confinement effect of the walls and behaves like a free particle. Thus, both ground and 
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excited states binding energies will become the free electron binding energies. 



IV. CONCLUSION 

In this study, we have carried out the deeply investigation of the electronic properties 
and the donor binding energies of a core/shell/ well/shell quantum dot structure. In order 
to determine the energy eigenvalues and corresponding wave functions, the Schrodinger 
equation is solved numerically with the shooting method in the effective mass approximation 
for a finite confining potential. The effect of a hydrogenic donor impurity located at the 
center of the multi-layered spherical quantum dot is also considered . The results for cases 
with and without a hydrogenic impurity are comparatively presented as a function of the 
core radius and the layer thickness. We have concluded that the attractive Coulomb effect 
of the impurity results in a significant change of the electronic structure of the multi layered 
quantum dot. The size of the barrier, well and core layers are also effective to determine 
the electronic properties of the system. Trends of all results are in a good agreement with 
other theoretical studies. We hope that this study will be very useful and contribute to 
understand the electronic properties of a core/shell/ well/shell spherical quantum dot with 
and without an impurity. 
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